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Abstract 

> 

I/-} ■ The basic mathematical assumptions for autonomous linear kinetic equations for a 

classical system are formulated, leading to the conclusion that if they are differential 
qq ■ equations on its phase space M, they are at most of the 2nd order. For open sys- 

\ terns interacting with a bath at canonical equilibrium they have a particular form 

\ of an equation of a generalized Fokker-Planck type. We show that it is possible to 

obtain them as Liouville equations of Hamiltonian dynamics on M with a particu- 
lar non-commutative differential structure, provided certain geometric in character, 
conditions are fulfilled. To this end, symplectic geometry on M is developped in 
• this context, and an outline of the required tensor analysis and differential geome- 

try is given. Certain questions for the possible mathematical interpretation of this 
structure are also discussed. 
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1 Introduction 



The understanding and description of irreversible evolution of macroscopic systems, pre- 
sumably towards equilibrium states is the aim of kinetic theory. This is done by formu- 
lating appropriate kinetic equations giving the time evolution of the state of the system 
which loosely speaking, is assumed to be some probability measure on the space of vari- 
ables describing the system, which henceforth will be called its phase-space. Observables 
are assumed to be well-behaved phase-space functions and experimental results refer to 
expectation values obtained via a bilinear (or sesquilinear) form on the cartesian prod- 
uct of the spaces of states and observables. Thus at the very heart of kinetic theory a 
probabilistic view point is rooted, the precise interpretation of which is often a matter of 
debate. 
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Theoretically speaking, kinetic equations are derived by following two different proce- 
dures: 



(i) stochastic arguments, based on some presumably plausible assumptions on the be- 
haviour of a large number of microscopic events characterising the system. 

(ii) application of more or less systematic approximation schemes on the exact micro- 
scopic dynamics of the system under consideration. 

Typical examples for (i) are the Fokker-Planck and Kramers equations (e.g. H ch.VIII). 
The first refers to the probability distribution function (p.d.f.) of the velocity of a heavy 
Brownian particle suspended in an equilibrium medium of light particles, the second one, 
to the p.d.f. of its phase-space position under the presence of an external field. The ba- 
sic assumptions for their derivation are that the microscopic dynamics of the Brownian 
particle are governed by Langevin's equation and that its p.d.f. is a Markovian diffusion 
process (see e.g. M ch.II, [[|). This procedure usually leads to model equations that 
are easier to study than the assumed more fundamental equations following from exact 
dynamics by (ii). Mathematically speaking, this procedure essentially replaces exact dy- 
namics by stochastic differential equations which in turn imply evolution equations for the 
p.d.f. of the system (e.g. Q ch.4, || ch.9). Though interesting and important from both 
the mathematical and physical point of view, it seems that one has to know somehow the 
kinetic equation he wants to derive and modify accordingly the corresponding dynamics. 
Moreover the modification is not easily interpreted physically (see the discussion following 
( ^.1[ ) below). That is, it seems that there is no general prescription of how this modifica- 
tion has to be carried out. In this paper, starting from mathematically totally different 
concepts we do arrive at results that may be interpreted stochastically, at the same time 
giving hints on the nature of such a general prescription (cf. the discussion following (|3.1|) 
below) . 

Typical examples for (ii) are the Boltzmann equation for dilute classical gases (or its 
quantum weak-coupling analogue, Pauli's equation) or the Landau and Balescu-Lenard 
equations for neutral plasmas, as they are derived from Liouville's (or von Neumann's) 
equations by using either iteration schemes and projection operator methods, or by using 
its equivalent, the BBGKY hierarchy of equations truncated on the basis of physical 
considerations.^ 

Although many specific equations can be derived by mathematically satisfactory (or 
even rigorous) methods in some particular limit of appropriate parameters of the system 
(for a survey see f33|), it is true that any approximation scheme, leading to satisfactory 



kinetic equations for particular classes of systems, runs into trouble as long as one tries 
to extend it to other systems and/or higher order of approximation (||12||, [|i~3] §5). For 



iSee e.g. § §§30, 41, @ ch.20, § ch.IX, § §2.4, §§4.3, 4.6, §VI. 
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instance the linearized Landau equation follows for spatially homogeneous plasmas as a 
2nd order approximation in the plasma parameter to the Liouville equation (plus some 
additional assumptions which need not be discussed here). Any effort to find its general- 
ization either in higher order of approximation and/or for inhomogeneous systems, runs 
into difficulties (e.g. equations violating the positivity of p.d.f or having no if-theorem, 
are obtained — see e.g. [H| in connection with [|TB[] §5), or involves highly ad hoc steps 
which are sometimes hidden in the calculations (see e.g. fl5| in connection with §2 below). 

Although one may be tempted to accept that there is no reason to expect that equa- 
tions in an approximate theory share all properties of the corresponding equations in 
the exact theory, in kinetic theory things are more complicated since the exact theory 
in this case is simply (classical or quantum) dynamics; however as it stands, microscopic 
dynamics of a system with a very large number of degrees of freedom is useless as a 
macroscopic theory since it does not incorporate irreversible evolution and its relation 
to a theory dealing with macroscopically defined quantities is remote, or at least not 
straightforward. Moreover from what has been said above, the situation is even worse, 
since it often happens that equations obtained at a lower level of approximation (e.g. with 
respect to expansion in some parameter) exhibit the correct properties, which however 
disappear in any higher level of approximation!^ 

In our opinion this is an inevitable consequence of the philosophy underlying kinetic 
theory, namely that irreversible evolution is simply an approximation to the exact (clas- 
sical or quantum) reversible dynamics. To put it differently, the aim of kinetic theory 
(description of macroscopic irreversible evolution) seems not to be compatible with its as- 
sumed fundamental laws (reversible microscopic dynamics), at least in a straightforward 
manner. Therefore it seems necessary that a fundamentally different approach to kinetic 
theory may not be worthless! 



2 Mathematical assumptions underlying kinetic the- 
ory and their implications. 

As already mentioned, at the very root of kinetic theory lies a probabilistic interpretation. 
Using this as a motivation and in order to clarify the difficulties noticed in the previous 
paragraph, we here describe the basic assumptions for a kinetic equation to be in prin- 
ciple, acceptable, and conditions under which an explicit general form can be obtained. 
We restrict the discussion to linear, autonomous kinetic equations, which include kinetic 
equations for the important class of open systems. To be more specific, we consider clas- 

2 See e.g. equations following from expansion of Liouville's equation or the so-called generalized master 
equation |(| §2.3, @ §2, fill §5, [^7) in connection with the last paragraph of the next section. See also 
for expansions of the Boltzmann equation. 
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sical systems, though similar results are known for quantum systems as well ( ||19|| , [ |16 
Theorem 4.2). 

We assume that 

(i) The phase-space M is a locally compact, Hausdorff topological space (e.g. the 
phase-space of a Hamiltonian system). 

(ii) Observables A are in C (M, C) the space of complex- valued continuous functions on 
M, having a finite limit at infinity. 

(iii) States £ are positive linear functionals on the observables, their values £{A) giving 
expectations. Since positivity of £, implies its boundedness in the supremum norm 



e.g. |21] p. 106-107), states belong to C*(M, C), the Banach dual with respect to 



this norm, which is the space of (regular) complex Borel measures on M. 

(iv) Kinetic equations for the observables have well-posed initial value problem, i.e. 
uniqueness and continuous dependence of solutions on the initial data hold. More- 
over, expectation values are continuous in time. These imply that solutions of a 
kinetic equation define strongly continuous semigroups of linear operators on the 
space of observables and that the corresponding adjoint equation defines such a 
semigroup on the state space.Q 

(v) The adjoint semigroup conserves positivity and normalization of the states, i.e. 
initial probability measures retain their character for all positive times. 

These plausible assumptions imply that the solutions of a kinetic equation for the 
observables define a Markov semigroup, i.e. a strongly continuous, positivity and normal- 
ization preserving one-parameter semigroup of operators on C(M, C) ||27|| . The termi- 
nology stems from the fact that such semigroups are in one to one correspondence with 
(time-homogeneous) stochastically continuous Markov processes described by a transition 
probability distribution p(t,x,E), which for each t,x is a regural probability Borel mea- 
sure on M (see e.g. [^3] p. 399, |24| Theorem 2.1, for outline of a proof). If we further 
assume that M is an n-dimensional differential manifold, the generator of the semigroup 
is defined on C 2 -functions and that a Lindenberg's type condition holds, 

lirn^ — — - = Xe(x) uniformly in x, 

where xe is the characteristic function of E, then it can be proved that the kinetic 
equation for observables (i.e. essentially the generator of the corresponding semigroup) 

3 Strictly speaking, strong continuity of the latter holds on a smaller subspace, which however uniquelly 
defines the adjoint semigroup Theorem 1.36). 
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has the formQ 

— = a ij (x)d i d j A + a\x)d i A , (2.1) 

where a % \ a 1 are continuous, a u is a non- negative definite matrix function and the sum- 
mation convention has been used, as it will be done in the sequel.^ Actually a % \a % are 
related to the diffusion and drift coefficients associated with the corresponding Markov 
process, given by the first two moments of p. 

In many approaches in kinetic theory which lead to satisfactory kinetic equations of 
the form (|2.1| ), the following condition holds 



a* = b* + Xjj , (2.2) 

where X l H is a Hamiltonian vector field. Moreover for open (classical or quantum) systems 
with Hamiltonian H, in interaction with "baths" in canonical equilibrium, the correspond- 
ing Hamiltonian function is an integral of the unperturbed motion of the open system, 
and therefore has in general the form —{H + F) with 

{F,H} = 0, (2.3) 

{ , } being the Poisson bracket, or the operator commutator, and the minus sign gives the 
correct Hamiltonian equation when the system does not interact with the bath. 
Moreover 

d ja ij - V = (3a ij djH , (2.4) 

where j3 is proportional to the inverse temperature of the bath.| These conditions will be 
used in section 4. 

Substitution of Q-Q in (|J gives 

BA 

— = -{H + F, A} + diicf'djA) - pa^djHfyA . (2.5) 

The essential conclusion drawn from the above discussion is (in a somewhat nonrigorous 
language) that linear autonomous kinetic equations for classical systems, which are differ- 
ential equations are necessarily at most of the 2nd order with nonnegative- definite leading 
coefficient and vanishing zeroth order term ( cf. eq. ( \2.1\ ) ) .[] 



4 See |25) Theorem 5.3 for outline of a proof, and |28] Theorem XIII. 53 for a partial generalization — 
the Levy-Khinchinc formula. 

5 For detailed proofs and more precise formulation of the various conditions see |27j] . 

6 See e.g. H Proposition 4.1, |§ eq.(2.19), § p.190, jl| eqs(III.26) and (111.16) together with 



(HI.19), m eq.(5.8). 



7 A corresponding result is also known for quantum systems pj|, 17 , pq. The special case of the 



Kramers-Moyal expansion of the linearized Boltzmann equation has been considered in 
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Since kinetic equations following from microscopic dynamics (see §1) are often differ- 
ential equations, the above result, severely restricts their form. In particular all methods 
based on power expansions of the solution of Liouville's equation with respect to some 
appropriate parameter, usually lead to unacceptable results since each approximation step 
increases the order of the differential operator by one.0 In the rest of this paper, we will 
present a different point of view, motivated by the discussion in the next section. 



3 Non- commutative geometry and stochastic calcu- 
lus. 

The usual approach to kinetic theory (method (ii) in Section 1) follows the scheme 

Microscopic (Hamiltonian dynamics) 

+ =>■ Kinetic equations 

Some systematic approximation scheme 

On the other hand, method (i) in Section 1, involves no systematic approximation scheme 
but a modification of microscopic dynamics which become now stochastic differential 
equations. In our opinion, the weak point in this case is that there is no general method 
of how to choose the microscopic stochastic differential equations. That is, it seems 
necessary that one somehow knows the kinetic equation he wants to derive and then 
writes down the corresponding stochastic equation. However it is not always clear how 
to interpret the latter. 

For instance it would be desirable to be able to derive by method (i) kinetic equations 
obtained from microscopic Hamiltonian dynamics. However, the latter in general involve 
operators with derivatives in the g-coordinates (see e.g. fl^| , [|13j, f28fl , p^). And it 




is a standard fact that stochastic differentil equations in phase-space involving Wiener 
processes X t of the form 

/ A q (q,p)dt + F q (q,p)-dX? \ 
" ^ A p (q,P) dt + F p (q,p) ■ dXf J ' 1 ' ' 

imply such derivatives in the corresponding kinetic equation if F q ^ (see e.g. || §4.3.3). 
However in many cases, one would like to interpret dq/dt simply as a velocity, all dynamics 
being incorporated in dp/dt (e.g. the velocity of a Brownian particle), in which case the 
above stochastic equation ( p.l|) is not easily interpreted physically. 

In the next sections we outline another approach in which the Hamiltonian character 
of microscopic dynamics is retained, but instead of approximation schemes, we make the 

8 See e.g. @ §2.4 eq.(2.199), § p.215, §IX.6 particularly p.280. 
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fundamental assumption that observables are now defined on a manifold with noncom- 
mutative geometrical structure. As it will be explained, this may be interpreted as a 
stochastic dynamical structure, though it is not known if this is necessary. Nevertheless, 
the formulation developped in the next section is closely related to what may be called a 
differential geometric approach to stochastic calculus^] though this relation will be further 
explored in another paper. 

The mathematical motivation for introducing noncommutative geometrical structure, 
stems from the fact that in the one-dimensional version of ( |3.1|) , X t being a Wiener process 
(i.e. Langevin's equation) the usual associative product between differential forms spanned 
by dt, dX t and functions of t, X t , implies that stochastic differentiation does not obey 
Leibniz's rule for the product of two functions, due to the appearance of 2nd derivatives in 
Ito's formula for the differential of such functions ( ||29|| §2, particularly eq.(2. 4)). However 
it is possible to define a modified noncommutative but still associative product between 
functions and differential forms such that Leibniz's rule holds ( [^9|| section 3). This prod- 
uct induces a noncommutative differential calculus on the ordinary algebra of functions of 
t, X t ( fl30fl section 2) via the basic commutation relations between functions and 1-forms 

[dt, t] = [dt, X t ]=0, 

[dX t ,X t ]=2 1 dt, { - } 

7 being the diffusion constant appearing in the usual Fokker-Planck equation obtained 
from Langevin's equation ( f2"9fl eq.(3.16)). Eq. (|3.2|) is readily generalized to a multi- 
dimensional Wiener process: 

[dt,t} = [dt,Xi} = 0, 

[dx\,xl}=Wdt, {6 - 6> 

where b % i is a symmetric bilinear form on the space of 1-forms.Q Therefore ( |3.3| ) can be 
rewritten more generally as 

[df,g] = b(df,dg)dt, (3.4) 

for functions f,g of X t ,t and a symmetric, bilinear form with components b^ in the 
"coordinates" X\ and such that dt lies in its kernel, i.e. b tt = b lt = 0. 

In view of the above discussion, we consider in the next section a (2n+ l)-dimensional 
manifold M x 1R and a differential calculus on a subalgebra A of the algebra of complex- 
valued functions onMxR satisfying Q3.4|) , and outline the formulation of extended Hamil- 
tonian dynamics as symplectic geometry on M x R. The basic result is that Liouville's 

9 See e.g. | p7[ , ps| particularly ch.VI, compare also with the approach in p9| ]. 

10 This is implied by the fact that since d satisfies Leibniz's rule, [dX t l ,X;?] = [dX|,JQ] which in fact 
shows that this commutator depends only on dJQ, dX| — see BG] §3. 
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equation for observables turns out to be of the form (\2.1\ ) with conditions $BJt ), ( \2.4\ ) hav- 
ing a simple geometrical meaning. Therefore it may be interpreted as a kinetic equation 
on the space of observables corresponding to a classical open system with phase-space M. 

The construction is coordinate-independent, and presupposes the definition of such 
fundamental concepts as vector fields, linear connections, symplectic structure and an- 
tisymmetric wedge product of forms on the differential calculus defined on A by ( J3.3Q , 
in close analogy with the corresponding concepts of the ordinary differential geometry.^] 
The derivation is formal in the sense that no systematic study of the representation theory 
of (|3.3| ) is made. Its already mentioned relation with stochastic calculus is a possibility, 
but it is not clear if others exist. A preliminary discussion of this problem is given in 
section 5. 

To make the presentation as transparent as possible, detailed calculations will be given 
in a subsequent paper in which tensor analysis for the corresponding noncommutative 
differential calculus is developped systematically. In fact this work constitutes only a 
first step towards a systematic formulation of kinetic theory as Hamiltonian (symplectic) 
dynamics in phase-space equipped with a noncommutative geometrical structure. 



4 Noncommutative symplectic geometry 

Let M be a 2n-dimensional manifold and A the algebra of smooth functions on M x M. 
The coordinate on IR will be denoted by t. Let Q be the universal differential envelope 
of A, i.e. fl is a ^-graded algebra Cl = 0, rg ^ f2 r with Cl r = {0} for r < and Cl° = A. 
Then there exists a linear mapping d : Q — > Cl of grade 1, which satisfies 

(i) dl = 0, where 1 is the constant function with value 1, 

(ii) d satisfies the graded Leibniz rule i.e. d{%pip') = (dip)ip' + (—l) r ip(dip'), for ip G fl r , 

(iii) d 2 = on all of Q and 

(iv) A and dA generate Cl. 

The universal differential envelope (Q,d) of A can be realized as follows (see Q): 
think of G Cl r as a function on (M x R)' r+1 , where for / G fi and a G M x R f(a) is 
the value of / as an element of A on a and for a , . . . , a r+ i G M x R and (ft G Cl r we set 

r+l 

(d(j))(a , . . . , a r+ i) := ^(-l) fc 0(a o , . . . , a fc _i, a k+1 , . . . , a r+ i) . 

fc=0 

11 For details on the systematic definition and presentation of general results in noncommutative geom- 
etry on a commutative algebra, see J3(l, particularly §§2,3. 
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Furthermore for G fl r , ip G fl s and do, . . . a r+s G M x 1R we set 

(#) («o, • • • , a r+s ) := [0(a o , . . . , a r )] [V>(a r , . . . , a r+s )] , 

for any non negative integers r, s. According to these rules 

(fdgh)(a,b) = f(a)[g(b)-g(a)]h(b) , (4.1) 

and hence df g 7^ gdf. On the ^4-bimodule of 1-forms fl 1 we define a new product 
• : fl 1 x fi 1 — > fi 1 as follows: for a, /? G fi 1 and a, 6 G M x R we set 

(am (3) (a, b) := a(a, b) j3(a, b) . 

Note that 

(fi a/2)* (giPg 2 ) = fi9i(a»(3)f 2 g 2 , 

and [df, g] = dfmdg. The universality of (ft, d) is expressed by the property that, if (fl, d) 
is any differential algebra on A, then there is a graded- algebra homomorphism tt : Q — > fl 
of grade such that tt\Q° = id^ and d o tt = tt o d (cf. [0 §3.1). 



Let b : fi 1 x fl 1 — > A be a symmetric left-right ^4-bilinear form i.e. b(f\af 2 , gi/3g 2 ) = 
figi b(a, 0) f 2 g 2 and assume that dt lies in the kernel of b. Let also X denote the differential 
ideal of fl generated by am (3 — dtb(a,/3), then we set fl := fl/X and it : fl —>■ fl for the 
canonical projection. Since X is differential the operator d : fl — > fl given by d = tt o d is 
well defined, and because dt lies in the kernel of b, a symmetric left-right ,4-bilinear form 
b : fl x fl — > .4. is uniquelly defined by 6 o (7r x 7r) = 6. 

Now set 

dfmdg:=[df,g], (4.2) 

and extend by left-right ^4-bilinearity, then it is easy to see that Tr(am/3) = a • (3, where 
a = tt(o) and similarly for (3. Obviously we have 

am(3 = dtb(a,(3) . (4.3) 

This is but a special case of the general procedure used to relate (fl, d) to any other 
differential calculus (fl, d) via an ^.-bimodule homomorphism tt, with X = ker tt and 
induce • on fl 1 by • on fl 1 (cf. |30] §3.2). Let £ l , % = 1, . . . , 2n be local coordinates on 



M then the elements of A can be written locally as functions of t, £ l , i — 1, . . . , In (see 
also section 5). If we set b^ := b(dC, 1 , d£ J ) we find the commutation relations — note that 
b(dt, dC) = b(dC, dt) = - 

[dC,e] = dtW , (4.4) 
[dt,t] = [dt,C} = [dC,t} = 0. (4.5) 
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These are special cases of 

[df,g] = dtb(df,dg). (4.6) 

By ((15|) we get for any f,g,heA that 

df • dg • dh = . (4.7) 

Applying d on a product of two functions and using ( [4.6|) we obtain 

d(fg) = (df)g + (dg)f - dt b(df, dg) . (4.8) 

Considering Q 1 as a right A module the dual module X is a left A module. We write 
(X, a) for the duality contraction. If we define 

Xf:=(X,df), (4.9) 



then we obtain from ( |4.8j ) 

X(fg) = g(Xf) + f(Xg) - b(df, dg)(Xt) . (4.10) 

The elements of X will be called vector fields. It can be proved that as a left A module 
X is free with basis given by d t , di, . . . , <9 2n , where 

8 =A 

1 ' ' (4.ii) 

4 := dt - . 

Thus for X G X we have 

X = X t d t + X i d i , (4.12) 

with X 1 := (X£ l ), X* := (Xt). More generally it can be proved that for any differential 
calculus on a differential manifold, satisfying (|4.7|) vector fields, i.e. elements of X, are 
second order differential operators without constant term, like ([4.121) . Thus the name 2nd 
order calculus is justified in this case (cf. section 5). As a further consequence fl 1 is free 
with dual basis dt, d£\ i — 1, . . . , 2n, and hence 

df = dtd t f + dCd l f. (4.13) 

Using the bilinear form b we define a linear mapping from Q 1 to X, a t— > a b by 

:=6(a,/9). 

Note that (dt) b = and (df) 6 = 
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Relations for forms of higher grade are obtained by applying d on equations 
3D. We find 

dW + d&d? = dtdbv , 

dtdt = , df dt + dtdC = . 1 ' ' 

These are special cases of 

dfdg + dgdf = dtdb(df 1 dg) , (4.15) 
which is obtained by application of d on (JO)). 



From (J75) and ( J^TW follows that all deviations of the present differential calculus 



from the classical differential calculus are proportional to dt. Therefore and by the second 
equation in t \4-14i ) it is also clear that for forms dttfi, for any <fi e Q all calculations proceed 
classically. This will help us to proceed more rapidly in what follows. 

We extend the • product to act between any 1-form a and an arbitrary form <p by 
using the "insert" opeartor _l of ordinary exterior calculus 

a m<f>:= dt a 1 'J0 . (4.16) 

On the right hand side everything is as in the ordinary differential calculus because of the 
presence of dt . It is not difficult to see that (cf. (|4.2|)) 



df.<P=[<PJ]. 

For a one form a the combination a* acts as a derivation of the product of differential 
forms, i.e. 

a • (0V) — ( a • 4 > )i ) + 4>( a • ^0 • 

The elements u of X which vanish on t, i.e. u(t) = are derivations of A and define a left 
A submodule X\ of X . With every u G X\ we associate mappings D u : Q — > Q defined up 
to terms lying in dt Q through the following relations 

D u (^) = (D u <P)^j + 0(A^) (mod dt) , (4.17) 

D u f:=u(f), (4.18) 

and 

D u dt = (mod dt) . (4.19) 

We write Di for Dq. and we set 

I),dr j -d^T',, (moddt), (4.20) 
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for the coefficients of the connection. For an r-form with 

1 , 



- ril ... v de i ---dr (moddt) 



we find 



where 



D k <f> = ^V fc n ... ir dC • • • dC - (mod dt) , (4.21) 



T ! 



V fc ir .. ir := d k 4> h ... ir - Y'uj.Oj..., Y 3 kir^ii-j ■ (4.22) 

Extending these definitions as usual to tensor products we obtain 

V i b jk := diV k + b^T^ef, . (4.23) 

In the following we demand the connection to be torsion free i.e. Y l [jk] — 0, from which 
follows 

d(f) = dC Di<f> (mod dt) , (4.24) 

and b- compatible, that is ViV k = 0. 

It should be emphasized here that in the context of the present noncommutative 
differential calculus it is possible to develop systematically tensor analysis so that the 
introduction of the above mentioned concepts of a connection and covariant derivative 
(cf. eqs( |4.17|) -( |4.20|) , (|4.21|) - (|4.23|) ) is made perfectly rigorous. However, this would lead 
us far away from our task to develop symplectic geometry and Hamiltonian dynamics and 
therefore it will be presented in a subsequent paper. Let us also remark here that in the 



old-fashioned index notation these differential geometric tools are first introduced in |41 
With the aid of these mappings and the • we define a new product in f2 

4>Ai[> :=<fnfi + (Di<f>)(d?*il>) . (4.25) 

It is easy to check that A is right A linear in both factors, i.e. 

(0/)A(^) = (0A^)^- (4.26) 

Furthermore it can be shown that the product is associative and as we shall see below 
also graded commutative, i.e. for e Q r and ip G £l s 

(Pa ip = (-irv A(p . 

An easy consequence of the above definition is 
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We define now an operator 



D/:=d/+idf«Ad/ 



(4.27) 



motivated by the fact that it satisfies the usual Leibniz rule 

B(fg) = (Df)g + (Dg)f 
Note that because of this property the 1-forms 



Df := df - -dtP 



(4.28) 



transform right-covariantly under a change of coordinates f J = f(£), i- e - we have D'f J ' = 
Df Here we have set P := V k r jk . Clearly dt, Df , . . . , Df n form a basis of ft 1 



with 



and 



Df = DCW) + dt(d t f) 



df = Df--dtV i V i d j f , 



where 



dtf ■= d t f - - 2 h"Y,0,f . 



(4.29) 



(4.30) 



The vector fields d t , <9i, . . . , dm form a basis of X dual to the above basis of ft 1 . From 
the definitions we find 



Df Df = df df + -dtdf F J 



(4.31) 



Df A Df = Df Df + dtdf b ej T l M 
and using ( |4.14j ) we obtain 



Df A Df + Df A Df = , 
dt A Df = dtdf , Df A dt = df dt , dt A dt = dtdt = . 



(4.32) 



(4.33) 



It is now easy to see using ( 4.26 ) and fl4.33|) , that A is antisymmetric. With the aid of 
the curvature of the connection P 



R % ikt '■— dk^ % n ~ deT\ k + P mfc r 



fcl j£ — U£L j k i" 1 m fcl j£ — 1 m£V. jk , 



m -pi ™ 

1 m/1 
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and the curvature 2-form 



n i j := -R l jM de^ e (mod dt) , 
we can prove the following usefull formulae 

DuDflf = -R k H] di\d k A(P) (mod dt) 



dDi 



-WiidjJcp) + d?T k ij D k <l) (mod dt) 



and 



d(df • 0) - dC •d(j) = -dt b ij D j( j) - dt j r jk (d£ k • 0) . 
Using these one can prove (see also [E^j §7) for (ft E Q r and t/> G ft 

d(0 Aip) = (d0) A V> + (-l) r A (dV>) + fi i i (6» i J0)(df • V) + dtU^D^D^) . (4.34) 
A special case of this formula is 

d(0/) = (d0)/ + (-!)> Ad/ + dt ^(D^djf) . (4.35) 



For a 1-form a = D£ l aj + dta 4 we find using ( 4.34 ) 



da = -Df A + dtd^ 



d t aci - diOi t - -b jk (V j V k a i + R e jki a e ) 



To connect the above results to symplectic geometry and Hamiltonian dynamics we need 
the exterior derivative of a 2-form u. In general 

uj = -DC A D£ j uu + dtdCuJi , 
2 



see (|4.33|) and the remark before ( |4.29|) ). Then a lengthy calculation gives 



du 



1 

3! 



Df A D£ j A 



-d[iUi k ] 



+ 



-dtdCdC 



Therefore a 2-form u is closed, if 



and 



3k] 



.(4.36) 



(4.37) 



d t ujij — d[iUij] — -b (VjtV^y — i? m H[i^j]m — R m kijUern) = . (4.38) 
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If we now assume as in ordinary symplectic mechanics, that dtUJij = then it can be 
proved that (|4.37|) , (|4.38|) imply that uJi satisfies 

dniufl + \b H V\ k Lo m ) = , (4.39) 

where only i,j are antisymmetrized. Therefore to takes the form 

to = ~Df A u tj + dtdC(diH - i& w V*wa) , 

for some function H . On the other hand ( |4.37| ) is the condition for u>ij to be closed in the 
ordinary exterior calculus, hence by Darboux's theorem it can be brought locally to the 
form 



I 



'.i 



Jij , J = (Jij) = [ .J Q I ' /= (<W), ot,/3= 1, . 



, = ~D£* A D? Jij + dtdddiH - ^-b M V k Ju) ■ (4.40) 



where J is the symplectic form in canonical coordinates. According to the interpretation 
of £ l in section 5 this can also be assumed in the present context. Therefore, to give our 
results a more familiar form we write u in this form 

However the subsequent calculations do not make any use of the fact that are taken 
here to be canonical coordinates. 

We already used the "insert" operator _l in the sense of the classical differential calcu- 
lus. It is necessary to extend its definition in the present context, since it is needed in the 
calculation of Hamiltonian vector fields. For X 6 X, f G A it is natural to put XAf = 
and as in the ordinary exterior calculus we set 

X-i<x-= (*,<*), UA1) 

X_l(0A^) := (X_I0) A^ + (-l) r 0A(X_l^) , V ' J 

for a 6 fi 1 , ^ G O r and if) £ fL If X = X % di + X t dt in the coordinate system in which u 
has the form ( |4.40| ), we have X_lD£* = X 1 and XAdt = X 1 . Therefore we find that 



X_iu = DC -JijX j + X\diH + Fi) + dtX^diH + F f ) , (4.42) 



Fi := -(l/2)^ fc V, J^ . (4.43) 



In order to write u in the form ( j4.40| ), it is necessary that it has maximal rank, which is 



2n since MxRis odd-dimensional. Therefore it has a 1-dimensional kernel given by the 
relation 

X_ito = . (4.44) 
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In ordinary extended Hamiltonian mechanics H is the Hamiltonian^ and the kernel is 
identified by definition with the space of Hamiltonian vector fields X. 
Equations ( |4.42| ), ( |4.44| ) give 

J i3 X j = X\diH + Fi) , X'idiH + Fi) = . (4.45) 

Setting J 1 - 7 := Jy we have J lk Jjk = <5j hence the first equation gives 

X i = -X'j^idjH + Fj) , 

and consequently the 2nd equation is identically satisfied. Therefore the Hamiltonian 
vector field defined by H is given by 



1 , 



X = X t ( (diH + F i )rd ] + (d t - -b tJ Vid. 



(4.46) 



As in ordinary extended Hamiltonian dynamics, the equation of motion for an observable 
A, i.e. for A G A takes the form X A = 0, which gives 



d t A = - \{H, A} + FiJ^djA] + -diib^djA) + -b ii T k ii d i A . 



(4.47) 



The Hamiltonian equation is now identical with the general kinetic equation (|2.5| ) for a 
classical open system in interaction with a large bath at canonical equilibrium, provided 
that 



a 



1 



b' J 



Ft 

ki 



d t F , 



-PdiH 



(4.48a) 

(4.48b) 
(4.48c) 



with {F, H} = 0. 

Condition flOgEQ is equivalent to d^F^ = which by flOS]) , flOSD , Q ) takes the 
form 



& (VfcV^Jij — R m kl[iJj]m — R m kijJlm) — • 



(4.49) 



When 6 U is nondegenerate, this is equivalent to the condition that the symplectic form 
Jij is harmonic with respect to the Laplace- Beltrami operator of b (see e.g. |55] p. 3). 
Equation (|4.48c| ) is equivalent to the condition that the canonical measure 



e = e-^d^ 1 ■ ■ ■ d£ 2n (mod dt) 



(4.50) 



12 Notice that by ( foi| ), ( fO^ ), eq.(^40|) is w = (l/2)dfd^J^ + dtdf strongly reminding con- 
ventional extended Hamiltonian dynamics. 
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is covariantly constant Die = (see e.g. |36j p. 215). 

The above results can now be summarized by saying that if the phase-space of a clas- 
sical open system E is endowed with a noncommutative geometrical structure ( |4.4j ), (f4.5|), 
because of its interaction with a bath at canonical equilibrium, then the corresponding 
Hamiltonian evolution of observables is identical to that given by conventional kinetic 
theory if the symplectic form is "harmonic" with respect to the "metric" connection de- 
fined by b % \ and the canonical (Maxwell-Boltzmann) measure defined on £ at the bath 
temperature is covariantly constant. Condition (|2.3|) is not expected to follow from the 
procedure followed so far, unless a precise relation of (|4.4j) , Qj.5|) to conventional dynamics 
is somehow made plausible. This will be considered in another paper. 



5 Discussion 

In the previous section we have given geometric conditions so that Hamiltonian dynamics 
in the context of noncommutative differential calculus defined on M x R by (|4.4j) - (|4.6|) 
can be interpreted physically as kinetic theory of classical open systems interacting with 
a large bath at canonical equilibrium. 

As already remarked at the end of section 3 the presentation so far is formal in the 
sense that the nature of the algebra A and its corresponding coordinate representation in 
terms of has not been specified. Here we discuss these questions further, but it should 
be emphasized that this is not done rigorously. Actually much remains to be done for the 
complete clarification of the problems addressed in this section. 

At the beginning of section 4 we remarked that a 1-form a in the universal differential 
envelope of the algebra of functions on a set N is a function a:iVxiV->C obtained by 
the obvious extension of (|4.1| ), that is of 

(f~dgh)(a,b) = f(a)[g(b) - g(a)]h(b) . (5.1) 

For the • we can show that Q5.1|) implies 

(d/i* • ■ • M/ r )(a, b) = (Mb) - A (a)) ■ ■ ■ (f r (b) - f r (a)) . (5.2) 

Relations like ai/3 — dt b(a, (3) = are in general incompatible with the above prescription 
of evaluating differential forms. Therefore if we do impose such relations (i.e. pass from 
(Cl, d, •) to (Q, d, •) as it is outlined at the beginning of section 4) and at the same time we 
still want to retain somehow an interpretation of •, similar to that given by (|5.2j ) then the 
elements of Q l cannot be functions on the whole of N x N. In fact such relations induce 
some structure on N x N by grouping together points of iV which may be considered as 
neighbouring. This is best illustrated by giving some examples. 

Take N := R the real line. Let x be the coordinate function on N and impose the 
relation dx • dx — dx = (cf. [0). If we want to keep the interpretation of one forms 
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as functions on some set Aq, this cannot be the whole of iV x N. It is obvious that Aq 
must be that subset of A^ x A^ on which the imposed condition is satisfied identically. If 
(a,b) E Ni G N x N then since x(a) = a we find 

= (dx • dx — da;) (a, b) = (b — a) (b — a — 1) . 

Hence in order for (a, b) to be an element of Aq either b = a or b = a + l. Hence 
Aq = {(a, a), (a, a + 1)| a G R}. This set gives a structure on the set iV by specifying 
which of its points are to be considered as neighbouring. Obviously the above condition 
specifies a discrete structure on R. The possibility to evaluate a one form on (a, b) with 
a and b not neighbours is still given, if b — a = m E and corresponds to the "integral" 

rb m 

/ a := ^2 a ( a + k — 1, a + k) . 

Ja k=i 

Applying the same reasoning to the relation df • dg = with smooth functions /, g : 
R — > C, then for a, b G R we find 

(df . dg)(a, b) = (f(b) - f(a)) (g(b) - g{a)) = f(x 1 )g\x 2 )(b - a) 2 = , 

where we have used the mean value theorem. Since this condition must hold for all smooth 
functions we must have e := b — a and e 2 = 0. Now this relation gives something trivial 
since e must identically vanish. But the relation df • dg = [df, g] = holds in the usual 
differential calculus and consequently it cannot be trivial. In fact one may interprete 
relation e 2 = as saying, that e is an infinitesimal of first order. In this sense R is again 
structured since now N% := {(a, a), (a, a + e)| a G R}. For arbitrary a < b G R the integral 
is defined by 

rb m 

/ a := lim^ a(x k -i,x k ) , (5.3) 
Ja k=i 

where this is obtained by taking the limit of vanishing width of the partition a = xq < 
xi < ■ ■ • < x m — b of [a, b]. Now the relation df»dg = integrated over [a, b] for arbitrary 
a < b G R gives 

rb m 

= / df • d# = lim J2 (/(a*) - f(xk-i)) (g(x k ) - #(x fe -i)) , (5.4) 
Ja k=i 

which can be expressed by saying that the "quadratic variation" of functions must van- 
ish. This is true if f, g are of bounded variation, a condition which is necessary for the 
Riemann-Stieltjes integral fa gdf to exist. 

For a second order calculus on smooth functions of one variable £, parametrizing 
A^ := R, we have by definition (section 4, eq.( |4.7| )) that df • dg • dh = 0, hence by 
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applying d£ • d£ • d£ = on (a, b) E R 2 we get (f (6) - £(a)) 3 = 0. Thus £(&) - £(a) is an 
infinitesimal of second order. Therefore N\ = {(a, a), (a,a + e), (a,a + e 2 ) \ a e £ _1 (R)}, in 
this case and consequently for given a e N we can move away from a in two ways, either 
to a + e or to a + e 2 . In this sense then N becomes structured and can be considered as 
2-dimensionalJ_j 

For a < b G R we define formally an integral as in ( |5.2p . Applying this on d£ • d£ • d£ = 
we obtain 

m 

lmX;(£(zfc)-£(z*-i)) 3 = (), (5.5) 
fc=i 



where again the limit is obtained as in (p]q). This relation can be expressed by saying 
that the "cubic variation" of £ must vanish. It is perhaps of independent mathematical 
interest to find equivalent characterizations of such functions. 

Since the quadratic variation of £ does not vanish in general, £ cannot be the usual 
coordinate function of R. Consequently we have R as a differentiable manifold but with 
a differential structure which is not the standard one. If we set dt := (l/6)d£ • d£ with 
some constant b and t a function on R, then by Q2.5|) dt • dt = and hence t is of 
bounded variation and can be taken to be the coordinate function on some copy of R. 
This additional coordinate t realizes somehow the fact that N is 2-dimensional. 

In the light of the above remark, if M — R 2n , it seems that should be interpreted as 
local coordinates defined on R 2n by an atlas IA not compatible with the usual one giving 
the standard differential structure of R 2n ; that is there are charts in IA not C fc -related to 
the identity mapping of R 2n , for some k > 0. Equivalently we may say that the identity 
mapping of R 2n , does not belong to U. More generally, for a 2n-dimensional differentiable 
manifold M the above discussion implies that if £ is a local chart of M in R 2n and £ a 
local chart of R 2n from U the functions / : M — > C belonging to A, are smooth functions 
of £ but are not smooth, not even of bounded variation, as functions of £, that is / o 
are smooth, but / o (£ o are not. 

The whole discussion in this section reminds us strongly of stochastic calculus on man- 
ifolds developped in the context of semimartingale theory, in particular when stochastic 
terms are given in terms of Wiener processes (see e.g. [[J7J, ||38|| ). In fact there are many 
results in stochastic calculus having an exact analogue in the formalism developped here. 
As examples compare (|4.8|) with (4) in |58 p. 134 and their properties, or elements of 



X eq. (|4.10|) with the characterization of 2nd order fields in [|38j Lemma 6.1. Moreover 
our relation of the "metric connection" with the drift term in the general Fokker-Planck 
type kinetic equation (f4.43|) suggests a close relation with the interpretation in stochastic 



calculus of a connection on a manifold as a mapping giving the "drift" of a 2nd order 



vector field (|37j p. 258-259). In fact it seems possible — and it will be examined elsewhere 
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The equality N = 1R is only set-theoretic. 



19 



- that our present model of noncommutative geometry can be realized in the context of 
stochastic calclulus. However, whether this is the only possibility remains an interesting, 
but to our knowledge, still unsolved problem. 
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